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ABSTRACT 

We present Hubble Space Telescope ACS and STIS FUV/NUV/optical imaging of the radio galaxy 3C 236, 
whose relic ~ 4 Mpc radio jet lobes and inner 2 kpc CSS radio source are evidence of multiple epochs of AGN 
activity. Consistent with previous results, our data confirm the presence of four bright knots of FUV emission 
in an arc along the edge of the inner circumnuclear dust disk in the galaxy's nucleus, as well as FUV emission 
cospatial with the nucleus itself. We interpret these to be sites of recent or ongoing star formation. We present 
photometry of these knots, as well as an estimate for the internal extinction in the source using line ratios from 
archival ground-based spectroscopy. We estimate the ages of the knots by comparing our extinction-corrected 
photometry with stellar population synthesis models. We find the four knots cospatial with the dusty disk to be 
young, of order ~ 10 7 yr old. The FUV emission in the nucleus, to which we do not expect scattered light from 
the AGN to contribute significantly, is likely due to an episode of star formation triggered ~ 10 9 yr ago. We 
argue that the young ~ 10 7 yr old knots stem from an episode of star formation that was roughly coeval with 
the event resulting in reignition of radio activity, creating the CSS source. The ~ 10 9 yr old stars in the nucleus 
may be associated with the previous epoch of radio activity that generated the 4 Mpc relic source, before being 
cut off by exhaustion or interruption. The ages of the knots, considered in the context of both the disturbed 
morphology of the nuclear dust and the double-double morphology of the "old" and "young" radio sources, 
present evidence for an AGN/starburst connection that is possibly episodic in nature. We suggest that the AGN 
fuel supply was interrupted for ~ 10 7 yr due to a minor merger event and has now been restored. The resultant 
non-steady flow of gas in the disk is likely responsible for both the new episode of infall-induced star formation 
and also the multiple epochs of radio activity. 

Subject headings: galaxies: active — galaxies: individual (3C 236) — galaxies: starburst — galaxies: jets 



1. INTRODUCTION 

Galaxies occupy a heavily bimodal distribution in color- 
magnitude space, wherein young, predominantly disk- 
dominated galaxies reside in a 'blue cloud' and evolve onto a 
chara cteristically quiescent, bulge-domin ated 'red sequence' 
(e.g.. iBell etail 12004 IFaber et alj 120071) . The underdensity 
of galaxies in the 'green valley' separating these populations 
implies that cloud-to-sequence evolution is swift, requiring 
a cessation of star formation more rapid than would be ex - 
pected in passively evolving systems (e.g. JCowie et al.ll996l) . 
Quasar- and radio-mode feedback models have been pro- 
posed as mechanisms by which star form ation may be trun- 
cated by the heating a nd expulsion of gas (ISilk & Reeslll998t 
Hopkins et ai] 120051: iCroton et al.1 120061: iSchawinski et all 



2006), as it is now known that quasar activity was two or- 
ders of magnitude more common at red shifts z ~ 2 than at the 
present time (e.g.. ISchmidt et ai][l9 91). This, considered in 
the context of de clining star formation rates in massive galax- 
ies at z ~ 2 (e.g. JPerez-Gonzalez et al.l l2008). along with the 
emerging consensus that most populations of galaxies har- 
bor quiescent black holes at th eir centers (hereafter BHs, e.g, 
Kormendy & Richstone 1995), has given rise to questions of 
whether all bri ght galaxies go throu g h one or more active 
phases (e.g., lHaehnelt & Reesl 119931; ICavaliere & PadovanH 



[19891) . In this scenario, the quenching of star formation via 
feedback from active galactic nuclei (AGN) ma y be one of 
the primary drivers of cosmic downsizing (e.g., ICowie etafl 
l!996t rScannapi eco et alJ l2005. and references therein). 

The relationship between the AGN duty cycle and the regu- 
lation of host galaxy stellar evolution is far more complicated, 
however, as it can play competing roles at successive stages 
of galactic evolution. AGN activity has been associated not 
only with quenching star formation on large scales, but also 
triggering it via ISM cloud compression from the propagat- 
ing relativistic jets associated with radio galaxies (e.g, the 
so-ca l led "alignment e f fect", [R ees 1989 1 iBaum & He ckman 
119891: lMcCarthvlll993t feest et aljboOOt iPrivon et alj|2008l) . 
Moreover, it is natural to expect a correlated (but not nec- 
essarily causal) relationship between AGN activity and star 
formation. The tight relationship b etween BH mass and 
host galaxy bulge velocity dispersion (iMagorrian et al1ll998h 
iFerrarese & Merrittll2000t iGebhardt et al]l2000l) implies that 
the growt h of the BH and the galaxy bulge are tightly cou- 
pled ( e.g jKauffmann & HaehnelfeOOOtlCiotfi & van Albadal 
1200 lb . It is therefore expected that, throughout the process of 
hierarchical galaxy formation, gas infall due to major merg- 
ers or tidal stripping from a gas-rich companion can fuel not 
only AGN, but also the growth of the hos t galaxy stellar 
component via infall-induced starbursts (e.g., di Matt eo et alj 
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FIG. 1. — The two radio sources associated with 3C 236. (a) 326 MHz WSRT radio contours (in green, from Mack et al. 1997) of the "relic" radio emission 
associated with 3C 236, overlayed on SDSS imaging of the same region of sky. The deprojected size of the jet is ~ 4.5 Mpc, making it the largest known 
radio galaxy and one of the largest objects in the universe, (b) Global VLBI 1.66 GHz radio con tours of the central 2 kpc CSS ("young") radio source from 
Schilizzi et al. 12001), overlayed on HST/SJIS NUV imaging of the star forming knots described by O'Dea et al. (2001). The jet axes of the Mpc- and kpc- scale 
radio sources are aligned on nearly the same position angle. 



120051) . A significant fraction (~ 30%) of nearby power- 
ful radio galaxies exhib it evidence of infall-induced star- 
bursts near their nuclei (|Smith & H eckman 1989; Alle n et al.l 
12002b iBaldi & Capettill2008b iTremblav et al.ll2009l) . suggest- 
ing that the phenomenon is both common and compera- 
ble t o the lifetime of the radio source (~ 10 7 — 10 8 yr, 
e.g., iParma et"aT1 [T999). The AGN/starburst connection 
is therefore likely real and fundamental to galaxy evolu- 
tion itself, and its characterization has been a major pur- 
suit of the past two decades (e.g. . iBarn es & H ernquistl 



19911 ISilk & Reeslll998b lFabianlll999b Id i Matteo et al.ll2005L 
Hopki ns et al.1 12005b ISoringel et al. | | 2005b ISilverman et al.1 
2008blOuillen & Bland-Hawthornl2008T) ~ 



A key discriminant in understanding the nature and evolu- 
tion of the AGN/starburst connection may be found in some 
radio galaxies whose morphology is clear evidence for multi- 
ple epochs of AGN activit y. Several such examples have been 
observed (e. g., 3C 219 -[Bridle et alfil^lClarke et al.ll992l 
0108+388 - iBaum et al.lll990l) . and have come to constitute 
a new class of "double-double" radio sources, representing 
~ 5- 10% of predominantly large ( > 1 Mpc) radio galax- 
ies (e.g., |Schoenm^ersetan|2Offi30). Double-doubles are 
characterized by outer ('older') and inner ('younger') radio 
sources propagating outwards amidst the relic of the previous 
epoch of activity. This apparently repetitive activity is thought 
to be a consequence of the AGN fuel supply having been in- 
terrupted, whether by exha ustion, smothering , or disturbance, 
at some time in the past dBaum et al.| [l990). This scenario 
is consistent with models of radio galaxy propagation (e.g., 
iKaiser et al.H2000t iBrocksopp et al.ll2007l) . The relative ages 
of the radio sources (and therefore the timescale over which 
the engine was cut off and re-ignited) can be estimated us- 
ing size estimates from radio maps coupled with a dynam- 
ical model for the jets and radio spec tral energy distribu- 
tions ( S EDs) of the radiating electrons (Schoenmake rs et al.1 
2000a b; O'Dea et al.ll2001l) . 



1.1. An important test case: 3C 236 



The nearby (z = 0. 1005) double-double radio galaxy 3C 236 
is an important test case in studies of the AGN/starburst con- 
nection, and is the basis of both this paper and a previous 
study by lO'Dea etall (120011) . 3C 236 is a powerful double- 
doubl e with a relic edge-brightened FR II (Fanaroff & Rilevl 
1974) radio source whose deprojected linear extent exceeds 4 
Mpc, making it the s econd largest known r adio galaxy (only 
J 1420-0545 is larger, Machals ki et alj|2008l), and even o ne of 
the largest objects in the universe dSch ilizzi etal]|200ll) . Its 
inner young Compact Steep Spectrum (CSS) source, whose 
apparent origin is cospatial with the nucleus, is only 2 kpc 
in extent and is morphologically reminiscent of a young clas- 
sical double. Anecdotally, 3C 236 was initially classified as 
a pure CSS source before it was associated years later with 
the massive relic FR II source (R. Laing, private communi- 
cation). The jet propagation axes of both the Mpc- and kpc- 
scale sources are aligned to within ~ 10° of one another (as 
projected on the sky). See Fig. Q] for radio contour over- 
lays of both sources, using 326 MHz Westerbork Synthesis 
Radio Telescope (WSRT) and 1.66 GHz Global Very Long 
Baseline Interfero metr y (VLBI) Network rad io mapping from 
iMack et al.1 (11997b and lSchilizzi et al.1 (12001b for the relic and 
CSS sources, respectively. 

In addition to its rare radio morphology, 3C 236 is also 
unique in that its nuclear dust complex is made up of an inner 
circumnuclear dusty disk that is somew hat misaligned with 
an apparently sepa r ate outer dust lane (|Martel et al.l 11999b 
IdeKoffetalJ 120001: ITremblav et all 12007b . The total dust 
mass in the complex is estimated to be ~ 10 7 M Q , based on 
Hubble Spac e Telescope (HST ) absorption maps and IRAS 
luminosities (de Ko ff et aLll2000b . In Fig.|2]we present a 1.6 
/zm / 0.7 //m absorption map of th e dust complex, originally 
presented in Tremblav et al. (2007) a nd made via divisio n of 
HST/NKMOS and W FPC2 data from lMartel etail (fl99l and 
iMadrid et ail (120061). respectively. 

The work bv lO'Dea et all d2001l) studied HST NUV and op- 
tical imaging of the central few arcsec of 3C 236, finding four 
knots of blue emission arranged in an arc along the dust lane 
in the galaxy's nucleus. Their original NUV data is presented 
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TABLE 1 

Summary of Observations of 3C 236 



Observatory Instrument Aperture Filter/Config. Waveband/Type Exp. Time [Orbits] Reference Obs. Date Comment 
(1) (2) (3) (4) ' (5) (6) (7) (8) (9) 



New Observations 



HST ACS HRC F330W {/-band Imaging 2516s [1] HST9S91 21 Oct 2003 SF Knots 

HST ACS HRC F555W V-band Imaging 2612s [1] HST 9897 22 Oct 2003 Dust Lanes 

HST ACS SBC F140LP FUV Imaging 6900s [3] HST 9897 21 Oct 2003 SF Knots 

HST STIS NUV-MAMA F25SRF2 NUV Imaging 2520s[l] HST 9897 19 Oct 2003 SF Knots 



Archival Observations 

HST WFPC2 PCl F702W ft-band Imaging 4 X 140s HST 5476 7 May 1995 Galaxy 

HST WFPC2 PCI F555W V-band Imaging 2x300s HST 6384 12Junl996 Galaxy & Dust 

HST STIS NUV-MAMA F25SRF2 NUV Imaging 1440s ffiST8275 03 Jan 1999 SF Knots 

HST NICMOS NIC2-FLX F160W NIR Imaging 1152s HST 10173 02 Nov 2004 Host Isophotes 




ARC SECONDS 



NOTE. — A summary of the new and archival observations used in our analysis. (1) Facility name; (2) instrument used for observation; (3) configuration 
of instrument used; (4) filter used; (5) corresponding waveband and specification of whether the observation was imaging or spectroscopy; (6) exposure time (if 
the observatory is HST, the corresponding number of orbits also appears in brackets); (7) corresponding reference for observation. If the observatory is HST, the 
STScI-assigned program number is listed; (8) date of observation; (9) comment specific to observation. 

the associated data reduction. In section 3 we present our 
results, including a comparison of our photometry with stel- 
lar population synthesis models, followi ng (in the interests o f 
consistency) the analysis strategy used in lO'Dea et al.l (120011) . 
We discuss our results in section 4, focusing on the role played 
by the AGN/starburst connection in the special test-case en- 
vironment of 3C 236. We summarize this work and provide 
some concluding remarks in section 5. Throughout this paper 
we use H = 71 km s" 1 Mpc" 1 , fl M = 0.27, and fl A = 0.73. 

2. OBSERVATIONS & DATA REDUCTION 

We present Hubble Space Telescope observations of the nu- 
cleus of 3C 236 obtained as part of the Cycle 12 GO program 
9897 by O'Dea and collaborators. These consist of U- and V- 
band imaging with the Advanced Camera for Surveys (ACS) 
High Resolution Channel (HRC), FUV imaging with the ACS 
Solar Blind Channel (SBC), and NUV imaging with the Space 
Telescope Imaging Spectrograph (STIS). Below we describe 
the specifics of each new observation presented in this paper, 
and provide a brief description of the archival imaging and 
spectroscopy that we also include in our analysis. We also de- 
scribe the steps taken to reduce the data. A summary of the 
new and archival HST observations utilized in this work can 
be found in in Table [TJ 

2. 1 . Cycle 12 ACS and STIS imaging 

Our HST/ACS observations were designed to enable high 
sensitivity multicolor photometry allowing for construction 
of an SEP of the b lue knots previously observed in 3C 236 
dO'Dea et al.ll200Tl) . The FUV, NUV, U-, and V-bands were 
chosen so as to provide multiple constraints on a young blue 
stellar population, while also enabling consistency checks and 
estimates on the amount of intrinsic reddening. At the redshift 
of 3C 236 (z = 0.1005), the U- and V- bands are located just 
blue-ward and red-ward of the 4000 A break, respectively. 
Exposure times were chosen to permit detection of the knots 
at adequate signal-to-noise (S/N) over a range of possible ages 
and intrinsic properties. Three orbits (6900 s) were obtained 
for the SBC FUV image using the F140LP long pass filter. 
We obtained one orbit each (~ 2500 s) for the HRC U- and 
V-band images using the F330W and F555W filters, respec- 
tively. The F140LP long pass filter on ACS SBC ranges from 



FIG. 2. — 1.6 /im 1 0.7 fjm colormap of the outer lane and inner dusty disk 
in the nuc leus of 3C 236 ma de vi a division of ' rt ST/NICMOS and WFPC2 
data from Martel et al. 1 1999) and Madrid et al. (2006), r especti vely. This 
absorption map was originally presented in Tremblay et al. 1 2007). 

in greyscale with black contours in Fig. [TJZ?). The lack of an 
obvious spatial relationship between the knots and the CSS 
source suggests t hat the starbursts ar e infall-induced rather 
than jet-induced. iMartel et alj (119991) had also detected the 
knots of emission in their HST /?-band imaging, albeit to a 
lesser degree as the knots are very blue. IP' Pea etalJ d200lh 
used their photometry in comparison with stellar population 
synthesis models to estimate upper limits to the ages of the in- 
dividually resolved star forming knots (seen in Fig.QJ>). They 
found disparate ages between the clumps of emission, find- 
ing two to be relatively young with ages of order ~ 10 7 yr, 
while the other two were estimated at ~ 10 8 - 10 9 yr old, com- 
parable to th e estimated age of the giant relic radio source. 
lO'Deaet all (120011) argued that 3C 236 is an "interrupted" ra- 
dio galaxy, and has undergone two starburst episodes approx- 
imately coeval with the two epochs of radio activity observed 
on Mpc- and kpc- scales. That work motivated follow-up ob- 
servations with HST at higher sensitivity and spatial resolu- 
tion, the results of which we present in this paper. 

We organize this work as follows. In section 2 we describe 
the new and archival data presented in this paper, as well as 
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FIG. 3. — (left) 2500 s V-band exposure of the nucleus of 3C 236, using HST/ ACS HRC with the broadband F555W filter. The outer dust lane is seen in 
white, while three of the four knots of star formation are seen to the south and west of the nucleus along the inner dusty disk, whose position angle is slightly 
offset from that of the outer.(rig/if) combined 6900 s HST/ACS SBC FUV (F140LP) image of the star forming knots observed along the inner dust structure of 
3C 236. The image has been smoothed with a two pixel Gaussian kernel. North is up, east is left. The two images are on the same scale. At a redshift of z ~ 0. 1, 
l"corresponds to ~ 1.8 kpc. 



~ 1350 A (with a hard cutoff) to ~ 2000 A with a pivot wave- 
length of 1527 A. The SBC Multianode MicroChannel Array 
(MAMA) has a spatial resolution of ~0"034 x 0."030 per 
pixel and a nominal field of view of 34"6x30."l. The SBC 
acheives a peak efficiency of 7.5% at 1250 A. The F330W and 
F555W filters on ACS HRC have central wavelengths (filter 
widths) of 3354 (588) A and 6318 (1442) A respectively. The 
HRC has a pixel scale of 0." 028 x0."025 per pixel and its field 
of view is 29"x26". It reaches a peak efficiency of 29% at 
- 6500 A. 

We have also obtained one orbit (2520 s) of NUV imaging 
with the Cs 2 Te MAMA detector on STIS. The F25SRF2 fil- 
ter has a central wavelength of 2320 A and a FWHM of 1010 
A, which permits geocoronal [o I] A 1302+ 1306 A contamina- 
tion in its bandpass, though its contribution is far lower than 
the detector background and is not expected to affect our re- 
sults. The F25SRF2 cutoff does not permit geocoronal Lya 
emission. The NUV-MAMA has a pixel scale of 0"024 and a 
field of view of 25" x 25". 

2.2. Archival data 

In this paper we make use of archival HST Wide-Field 
Planetary Camera 2 (WFPC2) imaging in V- and /?-band 
(F555W and F702W filters, respectively) obtained as part 
of t he 3CR snapshot p r ograms by Sparks and collabora- 
tors (IdeKoff etalJ [T996: McCarthy et all 1 19971; iMartel et al.1 



_S_|(_ 

[1999J). We also use the 1440 s STIS NUV-MAMA F25SRF2 
image of 3C 236, which formed the basis of the study by 
lO'Dea et al.1 (pOOlh. and was fo rmally presented as part of a 
data paper by lAllen et al.l (120021) . The //-band image obtained 
with the HST Near-Infrared Camera and Multiobject Spec- 
trograph (NICMOS2) in SNAP program 10173 (PI: Sparks 
iMadrid et al] 120061: iTremblav et all 120071: iFlovd et al.l 120081) 
was used for fits to the host galaxy isophotes. Imaging and 
spectrosco py from the Sloan Digital Sky Survey (SPSS) is 
also used (lYork et al.ll2000t lAdelman-McCarthy et al.ll2008l) . 

2.3. Data reduction 



The previously unpublished HST data presented in this pa- 
per were reduced using the standard On-the-Fly Recalibra- 
tion (OTFR) pipeline provided as part of the Multi-mission 
Archive at Space Telescope (MAST). For ACS, the OTFR 
pipeline combines calibrated and flatfielded dithered expo- 
sures using the multidrizzle routine with the default pa- 
rameters. The task calculates and subtracts a background sky 
value for each exposure, searches for additional bad pixels not 
already flagged in the data quality array, and drizzles the in- 
put exposures into outputs that are shifted and registered with 
respect to one another. From these drizzled exposures a me- 
dian image is created, which is then compared with original 
input images so as to reject cosmic rays on the drizzled me- 
dian. More info rmation on the specifics of multidrizzle 
can be found in Koeke moer et al.l (12002b . 

We have not performed post-pipeline processing on the two 
one-orbit ACS HRC drizzled images (U- and V-band), as 
for the purposes of this work we are concerned primarily 
with high surface brightness emission near the center of the 
galaxy, and the results from the OTFR pipeline were suffi- 
ciently free of cosmic rays and hot pixels to be deemed "sci- 
ence ready". For the three-orbit ACS SBC data, we combined 
the three individual calibrated and flatfielded files manually 
using mult idr i z z le with the default parameters. The sin- 
gle 6900s output image was left un-rotated with respect to 
North so as to avoid associated pixel interpolation errors in 
our photometry. Information regarding the reduction of the 
archival data utilized in this paper can be found in the appro- 
priate references cited in section 2.2. 

3. RESULTS 

In Fig. [3] we present the ACS/HRC V-band (at left in red) 
and ACS/SBC FUV (at right in blue) images in panels (a) and 
(b), respectively. The two images are aligned and on the same 
scale, with east left and north up. At the redshift of 3C 236, 
1 "corresponds to ~ 1.8 kpc. In Fig. Ufa) we present the same 
data in Fig. [3] as an overlay, with the FUV contours rendered 
in black on the V-band image for a clearer sense of the spa- 
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FIG. 4. — (a) ACS HRC V-band image of the nucleus of 3C 236 with ACS SBC FUV contours (sm oothed with a two pixel Gaussian). The knots of star 
formation along the inner dust structure have been labeled so as to be consistent with O'Dea et al. (2001). That work detected knots "3" and "4" as two distinct 
compact regions of emission in their STIS MAMA NUV ima ging, while we dete ct one continuous filament of emission in our more sensitive and higher spatial 
resolution FUV imaging. So as to maintain consistency with O'Dea et al. 12001) while making this clear we refer to this patch as "knot 3 + 4" throughout this 
paper, (b) The same HRC V-band image (with colors inverted to better emphasize the outer dust lane) with STIS S25SRF2 NUV contours. The knots of star 
formation are nearly morphologically identical in both the FUV and NUV, though the higher sensitivity of ACS is evident in the FUV contours at left, which 
highlight lower surface brightness features in the knots. Both panels are aligned and on the same scale. East is left, north is up. At the distance of 3C 236, 
l"corresponds to 1.8 kpc. 



tial relationship between the blue star forming knots (seen in 
the FUV) and the dust complex (seen in V-ban d). We have 
labeled each blue knot following the scheme in lO'Dea et al.l 
d2001l) to allow for easier comparis on of results. In the orig- 
inal STIS NUV image discussed in lO'Dea et all (1200 ll) (see 
Fig.QJ), knots "3" and "4" appeared to be separate, individu- 
ally resolved regions of emission. In our new higher sensitiv- 
ity and spatial resolution imaging with ACS we detect these 
two regions as one filament of emission that extends ~ 1". In 
the interests of consistency we have nonetheless named this 
region with the two labels "3" and "4" originally assigned in 
lO'Dea et al.l (120011) . and will hereafter refer to the filament as 
" Knot 3+4" when d iscussing both regions as a whole. As 
in lO'Dea et all d200Tl) . we will also refer to the "nuclear" 
FUV emission, seen in Fig. [4] as the contour cospatial with 
the nucleus in the underlying V-band. In Fig. |4jfe) we over- 
plot the STIS S25SRF2 NUV contours on the same V-band 
image (with colors inverted to highlight the continuum deficit 
due to the outer dust lane). The blue emission is morphologi- 
cally nearly identical in both the FUV and NUV, although the 
higher sensitivity (and longer exposure time) of the ACS data 
is evident in the FUV contours of Fig.HJa), which map lower 
surface brightness features than are seen in Fig.|4{Z?). We uti- 
lize the ACS/HRC [/-band image in our analysis, though do 
not present the image in a figure as it appears nearly identical 
to the knots in Figs. [3]^) and|4] 

We present an analysis of these new data in th e subsections 
below, framed in the context of past results from lO'Dea et al.l 
J200l . 

3.1. The outer lane and inner dusty disk 

In Fig. |3la) a significant deficit of galaxy continuum due 
to a dust lane extends ^10 kpc and runs in a northeast-to- 
southwest direction south of the nucleus at a position angle of 
~ 50°. A steep brightness and color gradient defines a ridge 
north of the outer lane, giving way to an inner circumnuclear 



disk of gas and dust whose major axis is oriented ~ 30° (off- 
set ~ 20° from the outer lane). It is not apparent whether the 
lane and disk are disjoint stru ctures or a continuo us, warped 
distribution of gas and dust. Ide Koff et al.l (120001) estimated 
the mass of dust in the complex to be ~ 10 7 M©, correspond- 
ing to a gas mass ~ 10 9 M w given the stand ard gas-to-dust 
ratio from the literature (Sod roski et al.|[l994l) . 

3.2. Properties of the star forming knots 

We have measured the total flux from each knot of star 
formation in the SBC FUV image (Fig. O?) as well as the 
STIS NUV (Fig. I!?) and HRC {/-band images. We have 
also measured the HRC V-band flux for knots 2, 4, and the 
blue emission associated with the nucleus. We have not per- 
formed V-band photometry for knots 1 and 3 as they are not 
clearly detected in this band (see Fig. Fluxes were mea- 
sured from the drizzled images using the apphot package in 
IRAF. The default drizzled pixel units for ACS are in electrons 
s , which scales to counts by a factor of the gain. As this 
is corrected for in the pipeline's calibration stage, "counts" 
and "electrons" can be considered equivalent for the purposes 
of this paper, regardless of the instrument being discussed. 
Count rates were summed in an aperture whose radius was 
chosen based on the size of the source being measured. Back- 
ground count rates measured through an aperture of the same 
radius were subtracted from the sum. Photometric conversion 
into flux units was applied by scaling the residual (source mi- 
nus background) value by the inverse sensitivity (the PHOT- 
FLAM keyword), converting the value from electrons s" 1 to 
flux units in erg cm" 2 s" 1 A -1 . Statistical ler uncertaintes were 
calculated from the measured count rates. The absolute pho- 
tometric calibration of the ACS SBC, ACS HRC, and STIS 
MAMA is 5%, 2%, and 5%, respective!^ Our photometry 



1 Maybhate, A., et al. 2010, "ACS Instrument Handbook" 
(Baltimore: STScI). 
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TABLE 2 

Photometry of the Blue Knots in 3C 236 



FUV Flux m F i40LP (FUV) m FS 25SRF2 (NUV) "1F330W ([/-band) 

Source erg s~' cnT 2 (mag) (mag) (mag) m F uv-m N uv >"FU V - »J(/-tand m NU v - m[/-band 

(1) (2) (3) (4) (5) (6) (7) (8) 



Corrected for Galactic Exctinction 



Knot 1 


8.305 X 10" 


18 


22.204 ± 0.025 


22.165 ±0.023 


22.074 ±0.015 


0.039 ±0.034 


0.130 ±0.029 


0.091 ±0.027 


Knot 2 


1.771 X 1Q- 


17 


21.382±0.017 


21. 840 ±0.019 


21.731 ±0.013 


-0.459 ±0.025 


-0.349 ±0.021 


0.109 ±0.023 


Knot 3 


1.134 X 10" 


17 


21. 866 ±0.021 


21.865 ±0.020 


21.498 ±0.012 


0.001 ±0.029 


0.368 ±0.024 


0.367 ± 0.023 


Knot 4 


1.670 x 10' 


17 


21.427 ±0.017 


21.788±0.019 


21.813±0.013 


-0.361 ±0.025 


-0.386 ±0.021 


-0.025 ± 0.023 


Knot 3+4 


2.836 x 10" 


17 


20.871 ±0.013 


20.961 ±0.013 


20.835 ± 0.009 


-0.090 ±0.018 


0.036 ±0.0158 


0.126±0.016 


Nucleus 


9.207 x 1Q- 


18 


22.092 ±0.024 


21.613 ±0.017 


20.359 ± 0.007 


0.480 ±0.029 


1.733 ±0.025 


1.254 ±0.018 


All 


9.037 x 10" 


17 


19.612 ±0.008 


19.703 ± 0.007 


19. 125 ±0.004 


-0.091 ±0.010 


0.487 ±0.008 


0.579 ±0.008 



Corrected for Galactic & Internal Extinction 



Knot 1 


4.240 x lCT 


16 


17.934 ±0.025 


18.141 ±0.023 


19.711 ±0.015 


-0.207 ± 0.034 


-1.777 ±0.029 


-1.570 ±0.027 


Knot 2 


9.044 x Kr 


16 


17.112i0.017 


17.817±0.019 


19.369 ±0.013 


-0.705 ±0.025 


-2.258 ±0.021 


-1.552 ±0.023 


Knot 3 


5.790 x 1Q- 


16 


17.596 ±0.021 


17.841 ±0.020 


19.136±0.012 


-0.246 ± 0.029 


-1.540 ±0.024 


-1.294 ±0.023 


Knot 4 


8.674 x 10" 


16 


17.157±0.017 


17.764 ±0.019 


19.451 ±0.013 


-0.607 ±0.025 


-2.294 ±0.021 


-1.687 ±0.023 


Knot 3+4 


1.448 x 10" 


16 


16.601 ±0.013 


16.937 ±0.013 


18.473 ±0.009 


-0.336 ±0.018 


-1.872±0.0158 


-1.535±0.016 


Nucleus 


4.701 x 10- 


16 


17.822 ±0.024 


17.589 ±0.017 


17.997 ± 0.007 


0.232 ±0.029 


-0.175 ±0.025 


-0.407 ±0.018 


All 


4.614 x 10" 


15 


15.342 ±0.008 


15.680 ±0.007 


16.762 ±0.004 


-0.337 ±0.010 


-1.420 ±0.008 


-1.082 ±0.008 



NOTE. — Photometry of the blue star forming regions in 3C 236. We have presented our results corrected only for galactic extinction, and then again with 
galactic and internal extinction corrections applied. The internal extinction correction was estimated from the Balmer decrement using SDSS spectroscopy 
(see section 3 and Table[3). lcr uncertainties have been derived from count rate statistics. The table does not include the systematic errors due to the absolute 
photometric calibration of the instruments, which are 5%, 2%, and 5% for the ACS SBC (FUV), ACS HRC ([/-band), and STIS MAMA (NUV), respectively. 
(1) Source name; (2) measured FUV flux of each blue source from the ACS SBC F140LP image; (3) apparent brightness of the FUV emission in the HST 
VEGAMAG system; (4) apparent brightness in the NUV; (5) apparent brightness in the [/-band; (6) FUV-NUV color; (7) FUV-[/-band color; (8) NUV-[/-band 
color. Larger (more positive) color values correspond to redder colors. 



w as corrected for galactic extinction using the scaling relation 
by lCardelli eUd](Tl989l) and a color excess E(B-V) = 0.011, 
as listed in the NASA/IPAC Extragalactic Database (NED). 
Summed and galactic extinction corrected counts were also 
converted to magnitudes normalized to the HST VEGAMAG 
system, defined such that the magnitude of Vega is zero in all 
Z/5Tbandpasses. 

As estimates on the intrinsic color of each knot is required 
in our analysis, and as each star forming knot is likely em- 
bedded in a great deal of dust, we have also corrected all pho- 
tometry for estimated internal extinction. We had originally 
proposed for and received STIS low dispersion long-slit spec- 
troscopy to allow measurement of the Balmer decrement (Ha 
/ H/3 flux ratio) in a 2"slit aligned with the knots, however a 
pointing error rendered the data unusable. I n lieu of this, we 
use th e measured Balmer decrement from Buttiglion e et al.l 
(2009), who presented line fluxes for 3C 236 using available 
Sloan Digital Sky Survey (SDSS) spectroscopy. 

Assuming an intrinsic Ha/H/3 theoretical line ratio of 
2.86 (i.e., a "case B" re combination scenario is assumed, 
lOsterbrock & Ferlandll2006l) . the observed Balmer decrement 
allows us to estimate the color excess associated with the in- 
te rnal extinction in the source, followin g the param eterization 
of lCardellietaLl (fl98l as described in Ate ketai] (120081) : 



2.5xlog(2.86/.R Qbs ) 
k{\ a )-k{\p) 



(1) 



where R \, s = F (Ha) / F (H/3) is the observed flux ratio, 
and the extinction curves at Ha and H/3 wavelengths are 
k(X a ) ~ 2.63 and k jXs) ~ 3.71, respectively, as given by 
ICardelli et alJdl989|). The Balmer decrement for 3C 236 mea- 
sured by Buttiglion e et al.l (120091) of /?„b s ~ 4.54 yields an 



E(B-V) Ha/H p w 0.465. The spatial resolution of the SDSS 
spectroscopy, while lower than that of HST, nonetheless al- 
lows for a reasonably believable estimate of the internal ex- 
tinction in the nucleus. In Fig. [5] we present the SDSS data 
for 3C 236, with a "zoom in" on the Ha and H/3 lines in 
Fig. EJfe). As for the purposes of this paper we are primar- 
ily interested in the spectroscopy solely for the Balmer decre- 
ment, we do not fit any of the l ines a nd instead use only the 
analysis from Buttigl ione et al.1 (120091) . We do note that a nar- 
row Ha absorption feature appears to be superimposed over 
the Ha+[Nll] lines, and possibly offset from the peak in the 
emission in velocity space. Ha absorp tion has been associ- 
ated with hot stars (Robi nson et al.|[l990l) . though we cannot 
confidently associate the feature with anything specific in the 
nucleus, as the spectrum is spatially blended over that size 
scale. 

Moreover, the lower spatial resolution of SDSS prohibits 
us from quantifying how patchy the extinction may or may 
not be on size scales of a few kpc (of order the scale over 
which the knots are distributed). It is not unreasonable to 
think that one knot may be more deeply embedded in the dust 
disk (and thus more significantly reddened) than another. Re- 
gardless, we expect that any variation in reddening between 
knots due to unresolved patchy extinction will be reasonably 
small when compared to the overall intrinsic reddening in the 
nucleus. Were this not the case, we would expect to see more 
severe color gradient shifts in the absorption map presented 
in Fig. |2] Instead, the colors of the outer lane and inner disk 
appear to be quite uniform over their projected lengths. 

In Table [2] we present the results from our photometric anal- 
ysis, including measured FUV flux, magnitudes, and colors 
for each knot. The data are presented first with only galactic 
extinction corrections applied, and then again with both galac- 
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TABLE 3 
Extinction Corrections 



Band 


Galactic Extinction 
A ( A )oal (mags) 


Internal Extinction 

A ( A )lta/H,a (mags) 


FUV (F140LP) 


0.1011 


4.270 


NUV (F25SRF2) 


0.0953 


4.025 


{/-band (F330W) 


0.0559 


2.362 


V-band (F555W) 


0.0292 


1.232 



NOTE. — Estimated and calculated corrections to photometry due 
to galactic and internal extinction . Galactic e x tinctio n estimated using 
E(B-V) = 0.011 and the law by Cardelli et al. 11989). Internal extinc- 
tio n calculated from the B aimer decrement F (Ha)/F (H/3) as measured 
by Buttiglione et al. 12009) from SDSS spectroscopy, using the method de- 
scribed in Osterbrock & Ferland (2006) and in §3.3. 

tic and internal extinction corrections applied. These correc- 
tions are given in Table 3. As the colors are always computed 
in magnitudes as "shorter wavelength"-"longer wavelength", 
la rger (more positive) values correspond to redder knots. As 
in lO'Dea"et al. (2001), we find that the knots vary relatively 
significantly in comparing colors that have not been corrected 
for internal extinction (columns 6, 7, and 8 in the upper half of 
Table 2). We note that knots 1 and 3 are significantly redder 
than 2 and 4. In correcting the colors using our estimate for 
internal extinction, however, we find that all knots are intrin- 
sically very blue, as expected. The sole exception is the UV- 
bright emission that is cospatial with the nucleus, which lies 
significantly redward of the "dust disk knots" in color space, 
even after the internal extinction correction has been applied. 
This could be due to an older stellar population or patchy ex- 
tinction leading to higher reddening along the line of sight 
toward the nuclear UV emission (or both). 

Again, we are unable to quantify whether or not there may 
be patchy extinction on scales finer than the SDSS spatial res- 
olution. As there is vastly more dust in the nucleus than in re- 
gions of the galaxy a few kpc outward from the dust complex 
(determined from our colormap in Fig.[2j, we find it likely that 
we have systematically undercorrected for internal extinction. 
However, a knot that is more heavily extincted than we have 
corrected for will be intrinsically bluer and therefore younger. 
Uncertainty therefore lies far more on the "young end" of age 
estimates than it does on the old, and we are abl e to estimate 
upper limits to the ages of the knots as done in lO'Dea et al.l 
( 120011) . We describe these results in the section below. 

3.3. Comparison of photometry with stellar population 
models 

We have compared our photometry with evolutionary mod- 
els from the stellar population synthesis code star burst99_] 
dLeifherer et al.lll999t I Vazquez & Leithererll2005l) . The pack- 
age incorporates Geneva evolutionary tracks as well as 
Padova asymptotic giant branch stellar models, and is now 
widely used for its effectiveness in accounting for all stellar 
phases contributing to the SED of a young stellar population 
from the FUV to the NIR. Several simulations were run, with 
parametric variations in stellar initial mass function (IMF), 
heavy element abundance, mass range, and whether or not the 
starburst continuously formed stars at a star formation rate 
(SFR), or whether it was instantaneous (wherein the starburst 
is modeled by a delta function, and the resulting population is 
allowed to age through time). 

2 http ://w w w. stsci . edu/science/starburst99/ 



Simulation results were ^-corrected to the redshift of 
3C 236 using its luminosity distance of ~ 457 Mpc. The cor- 
rected model spectra were then convolved through the rele- 
vant HST filter transmission curve using the photometric syn- 
thesis code SYNPHOlQin the STSDAS IRAF package. These 
corrected and convolved models could then be directly com- 
pared with our extinction-corrected photometry. 

An inherent challenge in comparing our data to predictive 
population synthesis models lies in the fact that we cannot 
know for certain what the absolute intrinsic color of each star 
forming knot is. As our internal extinction correction is based 
solely on the Balmer decrement as measured from (relatively) 
low spatial resolution SDSS spectroscopy, we are able to at 
best make a rough correction for internal reddening. As such, 
running several simulations with slight variations in model 
parameters was largely an exploratory exercise in which we 
qualitatively characterized how much a model might shift in 
color-color space given slight changes in, for example, the 
IMF or heavy element abundance. In general, these varia- 
tions were found to be small enough such that our uncertain 
photometry would be unable to discriminate between subtly 
different models. 

Ultimately, however, we are most interested in the rough 
ages of the knots as this study focuses on a possible coeval re- 
lationship between the blue star forming regions and the reig- 
nition of radio activity in the central engine. Model variations 
are such that our order-of-magnitude age estimates are largely 
independent of slight changes in model parameters. More- 
over, slight variations in our estimate for the internal extinc- 
tion would not affect the shape of the observed SED so much 
that we'd associate it with a significantly older or younger 
starburst. As such, we will focus our discussion on the "best 
fit" model we have found, with the caveat that all results pre- 
sented in this section are uncertain and highly dependent upon 
the intrinsic SEDs of the knots, which we can only roughly 
estimate. 

It is also important to note that there is line contamination 
in our broad-band filters that we have not corrected for, and 
it is necessary to quantify how much this might affect our re- 
sults, particularly for the FUV emission in the nucleus. In Ta- 
ble |4] we list intensities of selected bright emission lines (both 
directly measure d from SDSS spectroscopy as analyzed by 
Buttiglione et al.ll2009l and estimated using line ratios from 
Dopita et al.lfT997T measured in the M87 nucleus). In column 
(3) of the table we remark on whether or not the line falls 
within a passband used in our observations. Relative to this 
work, there is significant uncertainty associated with these 
line fluxes given the mismatch between SDSS resolution and 
fiber placement and the blue star forming regions we're inter- 
ested in. Moreover, there are even greater uncertainties in the 
estimated UV line strengths (as they are derived from ratios 
measured for the nucleus of M87, not 3C 236). We therefore 
use these fluxes only to roughly address the possible impact of 
line contamination on our results. There are relatively bright 
lines contaminating our optical and UV passbands, namely 
[Om]A5007 (in ACS HRC F555W) and CIVA1549 (in ACS 
SBC F140LP and STIS NUV-MAMA). However, correcting 
for this contamination would not change the overall shape 
of the SED relative to the comparison STARBURST99 model 
(which itself includes stellar and nebular emission, and as 
mentioned above, has been redshifted and convolved through 
the appropriate bandpass). Even correcting the photometry by 

3 http://www.stsci.edu/hst/HST_overview/documents/synphot/ 
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FIG. 5. — (left) Sloan Digital Sky Survey (SDSS) spectroscopy of the host galaxy of 3C 236. Significant emission and absorption lines have been labeled, and 
the positions of Ha and H/3, the lines most relevant for our purposes, have been marked with vertical dashed lines. The spectrum has been shifted to the rest frame 
using the known redshift of 3C 236 at z = 0.1. (right) "Zoom in" of the Ha and H/3 lines from the spectrum normalized to their local continuum. SDSS data are 
of suffi cient spectral resolution to resolve the two [Nil] lines at 6548 A and 6583 A about Ho at 6563 A (marked with a vertical dashed line). Buttislione et al. 
1 2009) measured the Balmer decrement = 4.54. See Table 4 for other emission line measurements. 



as much as a factor of two would not affect the age estimate 
we'd ultimately end up making, both for the knots in the dusty 
disk as well as for the nuclear FUV emission. We therefore 
do not correct our photometry for line contamination, given 
the uncertainties mentioned above, and because it ultimately 
does not affect our results. 

In Fig. [6] we present the main result of this paper, in which 
we compare the SEDs of the blue star forming knots as mea- 
sured in our data with the synthetic SEDs generated by STAR- 
BURST99. We plot four "snapshot" epochs of our best-fit 
model, nor malized t o a c ontinuous star formation rate of 1 
M yr" 1 , a S alpeteri (119551) IMF of slope a = 2.35, upper and 
lower mass limits of 100 M and 1 M Q , respectively, and so- 
lar abundances where Z = 0.020. These "snapshot" epochs are 
1 Myr (lowest line in blue), 10 Myr (second-to-lowest SED in 
green), 100 Myr (second highest in yellow) and 1 Gyr (top 
SED in red). As discussed previously, other star formation 
models were tested, but are not plotted in Fig. [6] as variation 
in overall SED slope and normalization was minimal. 

In black we overplot the measured SEDs of the blue star 
forming knots, as well as the blue emission in the nucleus, 
using our photometry as corrected for galactic and internal 
extinction. Individual knots are identified as per the legend in 
the upper right corner of Fig. [6] using the naming convention 
for the knots outlined in Fig. Hta). The black symbols on 
each SED mark the central wavelengths of the (from left to 
right) FUV, NUV, £/-, and V-band HST filters, respectively. 
As mentioned in the previous section, we were only able to 
obtain V-band flux measurements for knots 2, 4, and the blue 
emission associated with the nucleus (as they were the only 
regions clearly detected in the V-band image). 

3.3.1. The ages of the blue knots cospatial with the dusty disk 

We find the SEDs of knots 1, 2, 3, and 4 to best match 
the 10 7 yr model SED. Still, significant variation in color ex- 
ists among each knot, whose flux in a particular bandpass can 
vary by as much as a factor of five with respect to the others. 
Moreover, the SEDs compare somewhat inconsistently with 



the model (green 10 Myr line in Fig. |6j, as there is an appar- 
ent deficit of emission in the FUV, an excess in the NUV, and 
another deficit in the U-band (with respect to the model). 

These variations are due to (1) truly intrinsic differences be- 
tween the colors of the knots, (2) patchy extinction, and/or (3) 
globally under- or over-correcting for extinction. While we 
lack the ability to quantitatively discriminate between these 
scenarios, we find it naturally likely that a combination of all 
three work to produce the observed variations. It is also im- 
portant to note that we are studying UV emission in a very 
dusty environment that is optically thick to UV photons. The 
UV emission that we do see is therefore likely emitted on the 
near side of the clouds with respect to our line of sight, and 
we know little of the regions of the bursts that are more deeply 
embedded in the dusty disk. Nevertheless, we can confidently 
state that the observed SEDs of the knots are more consistent 
with a 10 7 — 10 s yr old starburst than they are with a 10 6 or 
10 9 yr old starburst (as modeled by STARBURST99, anyway). 

3.3.2. The age of the UV emission cospatial with the nucleus 

Along these lines, we find the nuclear emission to most 
closely match the 1 Gyr SED, though there is a deficit in its 
FUV flux with respect to the model (top-most red 1 Gyr line). 
Assuming that the FUV emission in the nucleus is indeed as- 
sociated with a 1 Gyr old starburst (an assumption we discuss 
below), it is possible (and even likely) that this FUV deficit 
is due to an undercorrection for internal extinction to the nu- 
cleus. Such an undercorrection would affect the FUV flux 
the most, while the V-band flux would remain relatively un- 
changed. The net result would be a steepening of the slope 
on the blue end, such that the SED of the nucleus might more 
closely correspond with that of the 1 Gyr model. Moreover, 
while we cannot quantify how patchy the extinction may be 
over these scales (as we have discussed several times previ- 
ously), it is reasonable to imagine the extinction being greater 
in the nucleus than it is on the edge of the dusty disk, given 
the inclination of the disk is more edge-on than face-on, such 
that the nucleus would be farther down the line of sight than 
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TABLE 4 

Emission Line Ratios & Radio Measurements 



Measurement Value 
(1) (2) 



In Passband? 
(3) 



Reference 
(4) 



Measured Intensities 


Log L(Ha) 


41.13 


No 




H/3/Ha 


0.22(4) 


F555W 




[Om]A5007/Ha 


0.57(2) 


F555W 




[Ol]A6364/Ha 


0.30(3) 


No 




[NlI]A6584/Ha 


0.69(1) 


No 




[Sll]A6716/Ha 


0.49(2) 


No 




[SII]A6731 /Ha 


0.35(3) 


No 




Log L(178MHz) 


33.56 


N/A 


2 


Log L(5 GHz) 


31.62 


N/A 


3 


Estimated UV Line Intensities 



Lya / Ha 


~6.7 


No 




ClV A1549 / Ha 


~ 1.2 


F140LP, NUV-MAMA 


4 


Hell A 1549 / Ha 


~0.3 


F140LP, NUV-MAMA 


4 


NIII]A1750/Ha 


~0.1 


F140LP. NUV-MAMA 


4 


Cm]A1909/Ha 


~0.3 


NUV-MAMA 


4 


Mgn A2798 / Ha 


~0.1 


F330W 


4 



REFERENCES. — (l)[B uttiglione et al. (2009); (2) Spinra d"eTaT] fl985l) ; 
(3) IButtiglione et aljf20iq) : (4) IDopita et alJ41997h 

NOTE. — (1) De-reddened logarithm of luminosity in erg s or line ratio 
(with respect to Ha); (2) measured value, with errors parenthetically pre- 
sented as percentages; (3) remark if the line (at the redshift of 3C 236) falls 
into one of the HST passbands used in this study; (4) reference. UV line 
intensities were estimated using measured ratios from the nucleus of M87. 
Significant uncertainties are associated with these estimates. 



the edge of the disk (Tremblay et all 120071) . Regardless, we 
do not expect that heavier extinction toward the nucleus (in 
comparison with the knots in the disk) would alter the shape 
of the nuclear SED so much that it would better correspond 
with a different age. 

Of course, it is important to estimate how much the AGN 
may contribute to the FUV emission in the nucleus. 3C 236 
is a low-excitation radio galaxy (LEG) that lacks broad lines 
(Buttiglione et al. 2009), suggesting that th e accretion region 
is ob scured along our line of sight (i.e., lUrrv & Padovanil 
1995). Contribution from the AGN to the nuclear FUV emis- 
sion would therefore come in the form of scattered light, 
which we do not expect to contribute significantly to the over- 
all UV excess (certainly not to a deg ree that wo u ld affe ct our 
order-of-magnitude age estimates). iHolt et ail (|2007|) simi- 
larly fit stellar population models to their ground-based SED 
for 3C 236, finding that a good fit required a red power-law 
component, suggestive that the UV excess is not likely associ- 
ated with scattered light from the AGN and is instead a young, 
reddened stellar population. Moreover, no point source is vis- 
ible in any HST imaging of the nucleus, and Spitzer Space 
Telescope IR spectroscopy obtained for 3C 236 {Spitzer pro- 
grams 20719 and 40453 by PI Baum and collaborators) does 
not show evidence for a hidden quasar continuum (although 
formal presentation of these data is forthcoming in future pa- 
pers). Hence, we expect the majority of the FUV emission in 
the nucleus to be due to stellar continuum from young stars, 
whose age we have estimated to be ~ 10 y yr old. 

3.3.3. The ages of the knots in context of other works 

To summarize, we conclude from Fig. [6] that the four star 
forming knots cospatial with the dusty disk are of order ~ 10 7 
yr old, while the stars in the nucleus are older, with an age of 



~ 10 9 yr old. In Table 5 we present the star formation rates 
associated with our "best fit" STARBURST99 continuous star 
formation model. We also list the mass range required for an 
instantaneous burst model to reproduce the observed fluxes, 
while noting that the continuous model is a better match to 
our data. Nevertheless, the mass ranges we estimate are con- 
sis tent with tho s e pred icted by the instantaneous burst models 
of lO'Dea et alj (12001b . Our estimated SFRs are also consis- 
tent with those derived from the O'Dea continuous star for- 
mation models , and are always on the order of a few M Q yr" 1 . 
See Table 3 in lO'Dea et al.l (1200 ll) to compare both their de- 
rived SFRs and mas ses with our results. 

Unlike this work. lO'Dea et al.l d2001 ) compared their pho- 
tometry to Bruzual & Chariot stellar population synthe- 
sis m odels dBruzual A. & Charlodfl99l ICharlot & Longhettil 
120011) absent an estimate for internal extinction (the required 
data was not available at the time). They found bimodality in 
the ages of the knots, wherein knots 1 and 3 were measured 
to be young, of order <; 10 7 yr old, and 2 and 4 were esti- 
mated at <J 10 8 - 10 9 yr old. In contrast, we estimate a nearly 
uniform age distribution among the knots in the dust disk of 
~ 10 7 yr old, while the FUV emission in the nucleus is likely 
to be ~ 10 9 yr. The differe nce in results between this work and 
that of lO'Dea et al.l ([200 1 ) likely arises from ( 1 ) different data 
sets (2) the lack of internal extinction correction in the O'Dea 
work (none was available at the time), and (3) comparison 
with Bruzual & Chariot models vs. STARBURST99 models in 
this work. Whatever the case, the age upper limits we esti- 
mate for o ur knots are consist ent with the "young" knots 1 
and 3 from lO'Dea et al.1 (120011) . That work also clearly em- 
phasizes that the older ages estimated for knots 2 and 4 are 
merely upper limits, and entirely dependent upon the intrinsic 
colors of the knots. 

The work by IHolt et all (120071) also estimated the ages of 
the star forming regions in 3C 236 by fitting young stellar 
population models to a ground-based SED from the ESO Very 
Large Telescope (V LT). Their results are consistent both with 
lO'Dea et al.l (120011) and this work, estimating similar ages for 
the stellar populations. In examining stellar absorption fea- 
tures, they further rule out models with ages > 1 .0 Gyr, as 
those models ov erpredict the depth of the Ca II feature. 

The work by iKoekemoer et al.l ( I 19991) performed a sim- 
ilar age dating exercise for the centrally dominant radio- 
loud elliptical galaxy in the cooling flow cluster Abell 2597. 
Much like 3C 236, that object possesses a compact radio 
source, a significant filamentary dust lane, and a network 
of clumpy knots and fi laments of blue continuum, which 
IKoekemoer et al.1 (119991) interpreted to be sites of recent star 
formation. Using single-burst Bruzual & Chariot models, that 
work derived ages for the knots of ~ 10 7 — 10 8 yr, comperable 
to the inferred age of the compact radio source. 

4. DISCUSSION 

4.1. Dynamics of the gas and dust in 3C 236 

As discussed in section 3.1 and seen in Fig. EJa), 3C 236 
possesses both an outer filamentary dust lane and an inner cir- 
cumnuclear disk who se semi-major axes are misaligned with 
one another by ~ 1 5° (iMartel et al.lll999tlde Koff et al.ll2000t 
lO'Dea et al.ll200l|) . These morphologies are strongly indica- 
tive of a galaxy that has recently acquired gas from a com- 
panion, suggesting that the outer lane is dynamically young. 
Below we motivate this assertion using simple dynamical ar- 
guments. 
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FIG. 6. — SEDs of the blue star forming knots as measured in our data, compared with synthetic SEDs as predicted by STARBURST99 stellar population 
synthesis models. We plot four "snapshot" epochs of one standard STARBURST99 model, normalized to a continuous star formation rate of 1 Mq yr~' , a Salpeter 
IMF of slope a = 2.35, upper and lower mass limits of 100 Mq and 1 Mq, respectively, and solar abundances where Z = 0.020. Other continuous star formation 
models were tested, but are not plotted here as variation in overall SED slope and normalization was minimal. In black we overplot the fluxes of the blue star 
forming knots, as well as the red nucleus, using our photometry as corrected for galactic and internal extinction (using the Balmer decrement). Individual knots 
are identified as per the legend in the upper right corner, using the naming convention for the knots described in 
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Gas and dust acquired through mergers or tidal stripping 
is ex pected to coalesce on a d ynamical timescale (~ 10 8 yr, 
e.g. lG^[l^7^lTubbl|[T980l) and precess about a symme- 
try plane of the host galaxy, finally settl ing into it on a pre- 
cession timescale of order a Gyr (e.g., iTohline et al.l 119821 : 
lHabe & Ikeuchi 119851) . During this time the gas will dis- 
sipate angular momentum and fall inward toward the nu- 
cleus at a rate dependent upon the structure of the poten- 
tial well and the star formati on efficiency of the gas (Barnes 
1996; Bekki & Shioval [l997l) . In the scenario proposed by 
Lauer et all d2005l) . filamentary distributions of dust that have 



not yet reached the nucleus would be classified as lanes, 
which might be thought of as transient structures that would 
eventually form a nu clear disk if given sufficient time (e.g., 
iTremblay et alJl2007l and references therein) . 

The recent study bv lTremblav et alj (|T007) lent evidence in 



support of this scenario. That work described a dichotomy 
between dusty lanes and disks in a sample of low-redshift 
(z < 0.3) 3CR radio galaxies (including 3C 236), finding 
round nuclear dusty disks to preferentially reside in round to 
boxy host galaxies, depending on their inclination with re- 
spect to the line of sight. Conversely, filamentary dust lanes 
which had not yet settled into the nucleus were found to reside 
exclusively in host galaxies with disky isophotes. Numerical 
simulations of dissipational mergers have shown that rotation- 
ally supported, disky systems are typically the result of gas- 
rich mergers, while boxy gala xies are often formed through 
dry (gas-poor) mergers (e.g.. | Barnesl [1996: Bekki & Shiova 
I1997L IKhochfa r & Burkerjl2005l) . Past studies of both radio- 
loud and radio-quiet ellipticals have shown that dust lanes 
are very often misaligned with the major axis of their host 
galaxy isophotes, while the opposite is true for dusty disks 
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FIG. 7. — (top left) //S77NIC2 1.6 (im image of 3C 236 with highlighted isophotal contours, (top right) 1.6 pm I 0.702 pm (H/R-band) colormap of the dusty 
disk in the nucleus of 3C 236, made via division of HST/NIC2 and WFPC2 data. Evidence for two disparate dust structures are seen in the colormap. The 
position angles of the outer dust lane and the inner disk-like structure are offset by ~ 10°. (bottom left) Semi-major axis position angles (P. A.) of ellipses fit 
to the //S77NIC2 //-band isophotes using the IRAF routine ellipse. The plot shows a several sigma twist in the isophotes between l"and 3", which is also 
visibly evident in the contours of the image at top left, (bottom right) Ellipticity of the fit isophotes from ellipse over the same scale. 

on the same scale). In Figs. Ul c ) an d (A) we plot isophote 
major axis position angle (RA.) and isophotal ellipticity, re- 
spectively, as a function of semi-major axis. These data 
stem from fits to the NICMOS //-band isop hotes using the 
IRAF task ellipse, o riginally performed bv|ponzelli et all 
(2007) and analyzed by ITremblav et all (120071) . Note from 
Fig- 13 c) that the isophotes are lopsided, particularly to the 
southeast between ~ 2"and 4"from the nucleus. The asym- 
metric isophotes are indicative of a stellar population that 
has yet to fully relax dynamically, suggesting a recent minor 
merger has taken place. 

Framed in the context of the results discussed above, the 
dual dust morphologies in 3C 236 present strong evidence 
of the host galaxy having recently (~ 10 9 yr) acquired gas 
from a companion. Following the scenario discussed by 
Trem blay et alJ (120071) . the inner dusty disk may have been 
native to the host galaxy before the recent gas acquision event, 



(e.g.. lTran et alJlTOQll and references therein). These results 
support a scenario in which nuclear dusty disks are native to 
the host galaxy pre-merger, while dust lanes are far younger 
structures, having recently been externally acqu ired through 
tidal stripping or a merger (ITremblav et alMool . 

3C 236 appears to contain both an outer dust lane and 
inner dusty disk. Moreover, the outer dust lane is slightly 
misalinged with the major axis of the host galaxy isophotes 
over a projected 4" linear extent centered about the nucleus. 
We mark the isophotal structure of the host in Fig. Iltfl), in 
which we have plotted isocontours on 1 .6 iim (analog o f H- 
band) KST/NICMOS imaging from iMadrid etail d2006h . In 
Fig. Ulb) we present a "zoomed in" view of the colormap 
originally presented in Fig. 12 with contours marking the mis- 
alinged outer lane and inner dusty disk. Note how the dusty 
disk appears to be better aligned with the major axis of the 
inner host galaxy isophotes than the lane (Figs.[7]2 and (b) are 
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TABLE 5 

Estimated Star Formation Parameters 



Source 
(1) 


Est. Required SFR 
(Mq yr _1 ,Cont. Model) 
(2) 


Est. Required Mass Range 
(Log Mq, Inst. Model) 

(3) 


Knot 1 


0.69 


7.3-8.7 


Knot 2 


1.47 


7.7-9.0 


Knot 3 


0.94 


7.5-8.8 


Knot 4 


1.41 


7.7-9.0 


Knot 3+4 


2.36 


7.9-9.2 


Nucleus 


0.47 


7.4-8.7 


All 


7.51 


8.4-9.7 



NOTE. — (1) Source identifier; (2) Estimated star formation rates (SFRs) 
using the "best fit" STARBURST99 "snapshot" epoch as determined based on 
Fig. |6]and as described in §4.2. We have chosen the 10 7 year SED as the 
"best fit" for all of the knots, and the 10' yr SED for the nucleus. The chosen 
STARBURST99 model for this estimate is normalized to a continuous star for- 
mation rate of 1 Mq yr~' , a Salpeter IMF of slope a = 2.35, upper and lower 
mass limits of 100 Mq and 1 Mq, respectively, and solar abundances where 
Z = 0.020. Other models were tested, but those results are not presented here 
as variation was minimal, and the SFR for each knot is generally of order 
~ 1 Mq yr -1 regardless of the model. (3) Required mass range of an instan- 
taneous starburst triggered 10 7 - 10 8 yr ago. The STARBURST99 parameters 
used are the same as in (2), except the starburst is modeled as a delta function 
at time zero, rather than forming continuously with a star formation rate. 

while the outer lane may be been externally acquired as a re- 
sult of that event, and is still in the process of migration toward 
the nucleus. 

As dust traces molecular gas accretion reservoirs that fuel 
AGN activity, this scenario works naturally wit h the notion 
that 3C 236 is an "interrupted" radio source. Ide Koff et alj 
(2000) estimated the mass of dust in the lane and disk to be 
~ 10 7 M Q , corresponding to a gas mass of ~ 10 9 M e . Such a 
significant gas mass could supply fuel to the AGN for a long 
period of time, allowing the radio galaxy to grow to its very 
large size (4 Mpc). The relic radio source also lends an unre- 
lated argument in favor of a scenario in which a large fraction 
of the gas mass was native to the host galaxy prior to what- 
ever event cut off or smothered the AGN activity for a period 
of time (ostensibly the same event responsible for depositing 
new gas into the system). Some significant fuel source was 
clearly necessary to power such a large radio jet, though it 
is not clear whether that original source is still present and 
observed in the form of the inner disk, or whether the origi- 
nal source was exhausted and the dust we observe now traces 
a newly deposited accretion reservoir. Imagining for a mo- 
ment that the inner dust disk is old and native to the host "pre- 
interruption", it is plausible to imagine the inner disk being 
perturbed, cutting off the BH fuel supply, then resettling back 
toward the nucleus after a short period of time (~ 10 7 yr), 
reigniting AGN activity. During this time the young, newly 
acquired dust from the "interruption" event would begin to 
settle into a filamentary lane, and begin migration inwards 
toward the nucleus on a much larger timescale (a few times 
~ 10 9 yr). 

The nature of this "interruption" event is not clear. The host 
galaxy of 3C 236 is in a very poor environment and has no ob- 
vious group members within 0.5 Mpc (b ased on a sta tistical 
correction for background contamination. IZirbell ll 997 ) . This, 
combined with the fact that the isophotes of 3C 236, while 
lopsided, are not highly irregular, suggests that the event was 
likely a minor merger. In this scenario, the presumably dust- 
rich and small donor galaxy would have been fully consumed 



by the much larger 3C 236. Were the event a major merger, 
one would expect to see more irregularity in the isophotes, as 
the "dynamically young" dust, the CSS radio source, and the 
fact that one can still detect the relic radio source even given 
synchrotron cooling timescales (~ 10 s - 10 9 yr) require the 
merger event to have happened relatively recently (no more 
than 10 7 - 10 8 yrs ago). A major merger of nearly equal- 
mass galaxies would likely result in structural irregularities 
that would last much longer than that timescale. 

4.2. Star formation in the nucleus in the context of the relic 
radio source 

Canonical estimates based on synchrotron cooling 
timescales would suggest that the relic ~ 4 Mpc FR II radio 
sourc e is not older tha n a few times 10 s yr (e.g.,[Parma et al. 
119991) . lO'Dea et afl (1200 ll) attempted to better constrain 
the age of the relic source, arguing that a lower limit to 
its age fmin.reiic could be estimated using simple dynamical 
arguments: 

fmi„, r eiic^7.8xl0 6 (^) ' yr (2) 

where vi„b e is the lobe propagation speed. Using a canonica l 
expansion s p eed of ~ 0.03c (e.g. jAlexander & Leahvll 19871) . 
lO'Dea et al.l d2001l) estimated the dynamical age of the relic 
source to be 2.6 x 10 8 yr. This would make for an old radio 
source, near to the time when synchrotron cooling would ren- 
der the lobes unobservable at higher frequencies as the elec- 
trons age. Of course, the relic source associated with 3C 236 
is among the largest in the universe, so the fact that it's likely 
old is not surprising. 

If the double-double radio source indeed arises from 
episodic activity in the AGN, at some stage in the life of 
3C 236 its nucleus would have entered a dormant phase and 
halted the collimation of its jets. This would have deprived the 
hot-spots of their energy supply on a timescale of order the 
nucleus-to-hotspot traversal time of the jet material ejected 
prior to the shut do wn of the AGN, assuming it advanced 
relativistically (e.g. , iBaumet alj Il990t iKaiser et al.l 120001; 
Schoenmakers et al. 2000a b). The production of "young 
electrons" in the hotspots by magnetohydrodynamical (MHD) 
turbulence or Fermi acceleration in the Mach disk is thought 
to cease when the hotspot is no longer fed by a jet, mean- 
ing the electron population will begin to a ge once the last of 
the remaining jet material has arrived (e.g.. lJones et al.lfl999L 
and references therein). If one believes radio spectral ageing 
techniques, then the ages of the youngest electrons in the lobe 
added to the jet traversal time from the nucleus to the hotspot 
should be of order the timescale o ver which the nucleus has 
been dormant. lO'Dea et al.l d2001l) estimated this timescale to 
be ~ 10 7 yr, an order of magnitude younger than the 2.6 x 10 s 
yr dynamical age of the relic source. 

Note that the ages estimated by these two techniques need 
not agree, as the ages of the youngest electrons in the lobe will 
correspond to the time when jet propagation ceased, but tell 
us nothing about how long the nucleus may have been active. 
Moreover, all of these age estimates in the radio are heavily 
dependent upon assumptions such as the propagation speed 
and the true source of young electrons. Indeed, the active 
phase corresponding to the creation of the relic radio source 
may have been far longer lived than the dynamical age of the 
relic source itself. We cannot know for certain how much of 
the lobes may have already cooled past the point of observ- 
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ability. Observations at lower frequencies may be enlighten- 
ing in this regard. 

As discussed in section 3.3.2, we have estimated the age 
of the nuclear starburst to be ~ 10 9 yr old, given the previ- 
ously discussed caveats and noting that it may be younger if 
we have undercorrected for internal extinction to the nucleus. 
Nevertheless, the approximate ages of the nuclear starburst 
and the ~ 4 Mpc relic source are just close enough to warrant 
noting that they might possibly be related to a common gas 
infall event. Having said that, we again stress that this is only 
a possibility, and we are unable to draw conclusions relating 
to such a connection in light of the significant uncertainties 
involved. 

4.3. The star forming knots in the context of the CSS radio 

source 

lO'Dea et al.l (120011) used arguments similar to those dis- 
cussed in the above section to estimate the age of the cen- 
tral CSS radio source, operating under the assumption that it 
was young and not "old and frustrated" (a reasonable assump- 
tion, especially in light of recent results on CSS sources, e.g. 
iHolt et aD l2009). Their estimate for its age ? m i n css was gi ven 
by 

f min , css ~3.2xl0 3 (^) ' yr, (3) 

where using the same lobe advance speed they found the age 
of the CSS source to be very young indeed, of order ~ l.Ox 
10 5 yr. As discussed in section 3.3, we have estimated that the 
star forming knots cospatial with the dusty disk are ~ 10 7 yr 
old. 

4.4. Is the recent episode of star formation coupled to the 

rebirth of radio activity in 3C 236? 

If the recent episode of star formation in the disk were di- 
rectly related to the event resulting in reignition of radio activ- 
ity giving birth to the CSS source, we would not necessarily 
expect a one-to-one correspondence among their ages. The 
dynamical time on scales of a few kpc, where the star form- 
ing knots are located, is far longer than the dynamical time 
on sub-pc scales where the AGN is fuelled. In this context, 
the difference in estimated ages of the star forming knots and 
the CSS source would constrain the timescale over which the 
gas is transported from kpc to sub-kpc scales. This can't be 
any longer than the age of the young stars (~ 10 7 yr), which 
corresponds not only with the dynamical time on kpc scales 
but also the estimated dor mancy perio d of the nucleus. 

The recent work by iWuvts et al.l (120 09) discusses the 
merger-driven models of Hop kins et al.l (120061) in the context 
of internal color structure in massive star forming galaxies. 
They find internal color gradients to be strongest during the 
merger phase, with blue star formation expected on scales 
outside a few kpc (approximately where we observe the star 
formation in 3C 236). In their previously mentioned work, 
iKoekemoer et al.l (119991) derived ages to star forming knots in 
Abell 2597 that were closely related to the inferred merger 
and AGN fuelling timescales giving rise to the compact radio 
source in that galaxy. 

Our observations have yielded similar results, and seem to 
suggest one possible "story" for the history of 3C 236. In con- 
sidering its double-double radio morphology in the context of 
its dynamically young dust complex and recently triggered 
compact starbursts, we suggest that 3C 236 has undergone 
multiple epochs of AGN activity due to a non-steady supply 



of fuel to the central engine. We suggest that the period of 
activity related to the 4 Mpc relic source was ended when the 
fuel supply to the central engine was cut off, whether by ex- 
haustion, strangulation, or disturbance. After a subsequent 
~ 10 7 yr dormant phase, infalling gas from a minor merger 
event reached kpc scales, where a starburst was triggered via 
cloud-cloud collisions amid collapse. The gas not involved in 
star formation reached the nucleus after a subsequent dynam- 
ical time, triggering reignition of the AGN and giving birth to 
the CSS radio source. 

We hesitate to make further suggestions relating to the star 
formation in the nucleus and its possible connection with the 
relic source, nor will we suggest whether or not the outer fil- 
amentary dust lane and inner disk are two distinct structures 
with different histories (i.e., recently acquired vs. native to the 
host pre-merger). We have argued in section 4.1 that discrim- 
inating between the two scenarios may shed light on whether 
the AGN "interruption" event was due to actual exhaustion of 
its fuel supply, or only dynamical disturbance of its fuel sup- 
ply. Whatever the case, the results of this work and those of 
lO'Dea et al.l (120011) strongly suggest that the transport of gas 
to the nucleus of 3C 236 has been significantly nonsteady over 
the past Gyr, giving rise to a unique galaxy that acts as an im- 
portant test case in studies of the AGN/starburst connection. 

5. SUMMARY AND CONCLUDING REMARKS 

We have presented follow-up HST ACS and STI S observa- 
tions o f the radio galaxy 3C 236, described by lO'Dea et alJ 
(120011) as an "interrupted" radio source. The galaxy is associ- 
ated with a massive relic ~ 4 Mpc FR II radio source (making 
it one of the largest objects in the universe), as well as an inner 
2 kpc CSS "young" radio source. This "double-double" radio 
morphology is evidence for multiple epochs of AGN activity, 
wherein the BH fuel supply is thought to have been exhausted 
or cut off at some time in the past, and has only recently been 
reignited. 

We present HST FUV, NUV, U-, and V-ban d imaging of 
four st ar forming knots, previously described bv lO'Dea et alJ 
(120011) . that are arranged in an arc along the outer edge of the 
galaxy's circumnuclear dust disk (which itself is surrounded 
by a misaligned outer filamentary dust lane). We have also 
detected blue emission cospatial with the nucleus itself. We 
describe these observations in detail, as well as the steps taken 
to reduce the data. We present photometry of the blue knots, 
and discuss our efforts to correct the data for internal extinc- 
tion to the source using the Balmer decrement available from 
archival SDSS spectroscopy. 

We compare the measured four-color SEDs of the star form- 
ing knots to synthetic SEDs from STARBURST99 stellar pop- 
ulation synthesis models, with the ultimate goal of roughly 
estimating the ages of the knots. We find that the four knots 
cospatial with the outer edge of the dusty disk are likely ~ 10 7 
yr old, while the FUV emission cospatial with the nucleus is 
likely older, at ~ 10 9 yr old (with the caveat that undercor- 
rection for internal extinction in the nucleus would lower this 
limit). We argue that the ages of the young knots are sug- 
gestive of a causal connection with the young central radio 
source. 

We frame these results in the context of 3C 236 as an ap- 
parently "interrupted" ra dio galaxy. Our resu lts are generally 
consistent with those of lO'Dea et al.l d2001l) . and we argue 
along similar lines that the transport of gas in the nucleus of 
3C 236 is nonsteady, wherein the active phase giving rise to 
the 4 Mpc relic source was cut off by exhaustion or distur- 
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bance of the AGN fuel supply. We suggest that this lead to a 
dormant period punctuated by a minor merger event, and the 
subsequent infalling gas triggered not only a new episode of 
star formation, but also ushered the galaxy into a new active 
phase giving rise to the young CSS radio source. 

Results such as these support the natural argument that in- 
falling gas largely fosters the relationship between the growth 
of the AGN and that of its host galaxy stellar component. 
In time, that relationship may evolve from mutual growth to 
the regulation of the latter by the former through quenching 
from AGN feedback. As infalling gas is a critical component 
in merger-driven hierarchical models, the triggering of AGN 
may be fundamental to galaxy evolution itself. In this way, the 
relevance of 3C 236 might extend beyond studies of episodic 
activity in radio galaxies to studies of AGN in the context of 
galaxy evolution as a whole. 
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